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1 Introduction. 



In this paper, we describe a new, systematic and explicit way of approximat- 
ing solutions of mixed hyperbolic systems with constant coefficients satisfy- 
ing a Uniform Lopatinski Condition via different Penalization approaches. 
In applied Mathematics like, for instance, in the study of fluids dynamics, 
the method of penalization is used to treat boundary conditions in the case 
of complex geometries. By replacing the boundary condition by a singular 
perturbation of the PDE extended to a larger domain, this method allows 
the construction of an approximate, often more easily computable, solution. 
We consider mixed boundary value problems for hyperbolic systems: 

d 

on {xd > 0}, with boundary conditions on {xd = 0}. The n x n real valued 
matrices Aj are assumed constant. Of course, we assume the coefficients to 
be constant as a first approach, aiming to generalize the results obtained here 
in future works. We assume that the boundary {xd = 0} is noncharacteristic, 
which means that detAd ^ 0. We denote by 
y := (xi, . . . , Xd-i) and x := Xd- The problem writes: 

Hu = f, {x> 0}, 

u\t<Q = , 

where the unknown u{t,x) G M", T : M" is linear and such that 

rg T = p; which implies that F can be viewed as a p x n real valued constant 
matrix. Let us fix T > once and for all for this paper. Let $7^ denotes 
the set [0,T] x M^J. and Tt denote the set [0,r] x W^'^. / is a function in 
H^{Q,^), g is a function in H^{Tt), where /c > 3 or /c = oo, such that: 
f\t<o = and g\t<o = 0. We make moreover the following Hyperbolicity 
assumption on Ti : 

Assumption 1.1. For all (r/,^) G M'^"^ x M — {0}, the eigenvalues of 

d-l 

Y^^jA^ + ^Ad 

are real, semi-simple and of constant multiplicity. 



(1.1) 
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Let us introduce now the frequency variable ( := (7, r, 77), where ir + 7, 
with 7 > 0, and r G M stands for the frequency variable dual to t and 
rj = (rji, . . . ,r]d-i) where rjj € M is the frequency variable dual to Xj. We 
note: 



Denote by M a x A^, complex valued, matrix; E_(M)[resp E-|-(M)] is the 
linear subspace generated by the generalized eigenvectors associated to the 
eigenvalues of M with negative [resp positive] real part. If F and G denote 
two linear subspaces of C-^ such that dimF + dimG = N, det{¥, G) denotes 
the determinant obtained by taking orthonormal bases in each space. Up to 
the sign, the result is independent of the choice of the bases. We shall now 
explicit the Uniform Lopatinski Condition assumption: 

Assumption 1.2. (H, F) satisfies the Uniform Lopatinski Condition i.e for 
all C, such that 7 > 0, there holds: 



The mixed hyperbolic system (II. ip has a unique solution in 
H^{Q^), and, since Ti is hyperbolic with constant multiplicity, for all 7 
positive, the eigenvalues of A stay away from the imaginary axis. More- 
over, as emphasized for instance by Chazarain and Piriou in [3j and Me- 
tivier in [8], there is a continuous extension of the linear subspace E_(y4) 
to {7 = 0, (r, r/) 7^ Ojjd} that we will denote by E_(A). E+(^) extends as 
well continuously to {7 = 0, (r, r/) 7^ O^d} and we will denote E_|_(^) this 
extension. Moreover, there holds: 



We can refer the reader to [3], [6], [7j, or [8] for detailed estimates concern- 
ing mixed hyperbolic problems satisfying a Uniform Lopatinski Condition. 
Moreover, we can refer to ^lOj for the proof of the continuous extension of the 
linear subspaces mentioned above in the hyperbolic-parabolic framework. 

Remark 1.3. As a consequence of the uniform Lopatinski condition, there 
holds, for all C : 




(1.2) 



|det(E_(A),kerr)| > C > 0. 



E^(^)0E+(A) = C- 



rg T = p = dimE_(y4((^)). 
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1.1 A Kreiss Symmetrizer Approach. 



We will now describe a penalization method involving a Kreiss Symmetrizer 
and a matrix constructed by Ranch in [12j , in the construction of our singular 
perturbation. Note well that we have some freedom in both the choice of the 
Kreiss Symmetrizer and of Ranch's matrix. Let us denote respectively by u, 
f, and g the tangential Fourier-Laplace transform of u, f, and g. Since the 
Uniform Lopatinski Condition is holding for the mixed hyperbolic system 
(jl.ip . there is, see [9J a Kreiss symmetrizer S for the problem: 



That is to say there exists a matrix S{(^), homogeneous of order zero in 
C~ in M+ X M'^ - {Ojjd+i} and there aie X > 0, 5 > and Ci such that: 

• S" is hermitian symmetric. 



An algebraic result proved by Ranch in |12| can be reformulated as follow, 
and a proof is recalled in section 12. 2t 

Lemma 1.4. There is a hermitian symmetric, uniformly definite positive, 
N X N matrix B such that: 



Moreover B depends smoothly of C,. 

Remark 1.5. This result is proved by constructing explicit matrices satisfy- 
ing the desired properties. Thus, it is not merely an existence result and we 
can use the explicitly known matrix B in our construction of a penalization 
operator. 

Let us denote by R := B^ and Sr := R-^SR'^. We will denote by 
the projector on E_(S'k) parallel to E_|_(5'ij) and by P+ the projector on 
E+(S'j?) parallel to E_(5r); and P"*" denoting the associated Fourier 
multiplier. We recall that, denoting by the tangential Fourier transform, 
the Fourier multiplier P~(5j, dy,^) [resp ^'^{dt,dy, 7)] is then defined, for all 
w G H''{R'^+^), and 7 > 0, by: 



(1.3) 




{x > 0} 



• K (SA) > Xld. 



• S>5Id- CiT*T. 



kerF = E+((5)-i5). 
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[resp 

in the future we will rather write: 

We fix, once and for all, 7 > big enough. Let us consider then the solution 
of the well-posed Cauchy problem on the whole space (jl.4p : 

Ult<o = 0, 

where 

M := -e'^^AdSr^RFTRe^'^^ 

and S_{dt,dy) [resp R{dt,dy)] denotes the Fourier multiplier associated to 
5(C) [resp R{C)]- Let us define g by: 

9-=e g. 

In what follows, g will denote the Fourier-Laplace transform of g. Let us 
denote by 

U := vTlxKO + ulj;>o = u"lx.<o + ulx>o- 

u denotes the solution of (jl.ip . and thus belongs to H''{^^). u~ is a function 
belonging to H^{Q.7p) and such that u~|x=o = u\x=o- More precisely, u~ can 
be computed by: e'^^T~^ {R~^{v~ + P~F^)) , where v~ is the solution of the 
problem: 



r SRdxvr - F+SrArv- = F+SRARF-Tg, {x < 0}, 
\ v~\x=o = F^Ru\x=o, 

and u denotes the Fourier-Laplace transform of the solution u of (jl.ip . 

Theorem 1.6. For all k > 3, if f e i7''(J7j) and g G H''{Tt), then there 
holds: 

where vf denotes the solution of the Cauchy problem (|1.4p and n denotes 
the solution of the mixed hyperbolic problem (II. ip . If g = then: 

\\u'^ — u~\\ ,, 3 + — n|| ,, 3 , = 0(e). 
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Of course, since is defined for all {x E M}, its limit as e ^ 0"*", u is 
can be viewed as an "extension" of u on the fictive domain {x < 0}. The 
"extension" resulting from our method of penalization gives a continuous u 
across {x = 0}, while the method used in [2] gave simply: :u|a;<o = 0- We 
have the following Corollaries: 

Corollary 1.7. Assume for example that f € and 
g G H°°{Tt) then 

M - u\\Hs(n+) = ^(e); > 0- 

Corollary 1.8. /// belongs to L^(J7y) and g = then: 

lim - u\\l2,^ ) = 0. 

One of the interest of this first approach lies in the rate of convergence 
of vf towards u. Indeed, in general, a boundary layer will form near the 
boundary in this kind of singular perturbation problem. For example in the 
paper by Bardos and Ranch [2], as confirmed by Droniou 4j, a boundary 
layer forms. It is also the case in [TT], as analyzed in our Appendix. There 
are also boundary layers phenomena in the parabolic context: see the ap- 
proach proposed by Angot, Bruneau and Fabrie [T] for instance. However, 
surprisingly, and like in the penalization method proposed by Fornet and 
Cues in [5], our method allows the convergence to occur without formation 
of any boundary layer on the boundary. As a result, this leads to the kind 
of sharp stability estimate given in Theorem II. 6i These results concern the 
case where / and g are sufficiently regular. The reason is that we construct 
an approximate solution. In the case of g only in L^(Tj'), such a simple 
treatment does not work. However, let 5 > be given. If we approximate 
/ and g by smooth functions G H°°{i}^) and g^, G H°°{Tt) such that 
11/ ~ /!^IIl2(q+) < S and \\g — 5i/||l2(Xj,) < S, by the uniform Lopatinski 
condition, we get: 

where is the solution of the mixed hyperbolic problem (jl.ip with data 
and gu- We can now apply Corollary 11.81 to u^, and obtain by penalization 
a sequence in L^(Qt) such that: limj^o+ '^l = '^i^ L'^{Q^). Finally, by 
choosing, e sufficiently small, we get \\u — uI\\^2(^q+^ < 2C6. By choosing e 

and ly as functions of 5, and noting u^^) = "U^^^l , we have: 
(1.5) \\u^s) - u\\L^in+) = 0- 
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1.2 A second Approach. 

In the first approach we have just introduced, it is necessary to compute 
a Kreiss's Symmetrizer and a Rauch's matrix. In view of future numerical 
apphcations, we wih now introduce another method preventing the compu- 
tation of these matrices. The price to pay is that we need the prehminary 
computation of v, which is by definition the solution of the Cauchy problem 
on the free space: 



(1.6) 



'Hv = /, {t,y,x) G Qt, 
v\t<o = V(y,x)GE^. 



Let us denote P (C) the spectral projector on IE_(A(C)) parallel to 
E+(A(C)), and P+(C) the spectral projector on E+{A{C)) parallel to E_(^(C)). 
Let us introduce P^(9t, Sj^, 7), the Fourier multiplier associated to P^(C). 
Let us denote by 11 the projector on E_(j4((^)) parallel to KerF, which has 
a sense because of the Uniform Lopatinski Condition and denote II the 
associated Fourier multiplier. We define then h by: 

h := e-' (P-(e-^*7;|,=o) + ne-^*(g - v\,=^)) , 

where g denotes the function involved in the boundary condition of the 
mixed hyperbolic problem (jl.ip . Now, let us consider the following singularly 
perturbed Cauchy problem on the whole space: 

[ nu' + -Ade^'p-e-^'un^^o = /lx>o + -vlde^*/il,<o, 
(1.7) I e e 

[u'\t<o = . 
Let us denote by 

u := u~lx<o + ulx>o = u~lx<o + ulx>o- 

u denotes the solution of (jl.ip thus belonging to H''{Q,^) and u~ is a function 
belonging to H^{Q.^) and such that u~\x=q = u\x=o- More precisely, u~ can 
be computed by: e^*T^'^{T{h)+v~'), where v~ is the solution of the problem: 

= 0, {x< 0}, 




and u denotes the Fourier-Laplace transform of the solution u of (jl.ip . The 
problem (jl.7p is well-posed and, for all e > 0, there exists a unique 

g H^{Q,t) solution. We will fix 7 adequately big beforehand. We observe 
then the following result: 
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Theorem 1.9. For all k > 3, if f e H''{n:^) and g £ H^{Tt), then there 
holds: 

\W - ll//fc-3{n-) + 11^'' ~ ^ll//'=-3{n+) = ^(^)' 

where u'^ denotes the solution of the Cauchy problem (|1.7p and u denotes 
the solution of the mixed hyperbolic problem (jl.ip . 

The singular perturbation involved in the definition of u'^ does not de- 
pend either of Kreiss's Symmetrizer or Ranch's matrix. As a result, for 
this method of penalization far less computations are necessary in order to 
obtain our singular perturbation. Note well that the proof of the energy 
estimates in Theorem 11.91 is completely different from the proof of the en- 
ergy estimates in Theorem 11.61 Indeed, for our first approach our singularly 
perturbed problem was treated as a Cauchy problem, contrary to our second 
approach where it was interpreted as a transmission problem. 

Corollary 1.10. Assume for example that f E H'^{il^) and 
g E H°°{Tt) then 

Of course, we see that the same problem of regularity arises in Theorem 
11.91 and Theorem 11.61 However, by a simple density argument, we can also 
prove here the exact analogous of (jl.Sp . 

Remark 1.11. In the case where / = 0, then the solution v of ()1.6p is 
V = Q and thus, the perturbed cauchy problem (ll.7p rewrites: 

' Hu' + - Arfe^*p-e-^*n^ U<o = -A^e^'e-^' {ike-^'g) U<o, {x 
< e e ' 

^ u^\t<o = . 

2 Underlying approach leading to the proof of 
Theorem II. 6L 

2.1 Some preliminaries. 

Since the Uniform Lopatinski Condition holds, there is S, homogeneous of 
order zero in and such that there are A > 0, 5 > and Ci and there holds: 

• 5 is hermitian symmetric. 

• ^ {SA) > Xld. 
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• S>6Id- CiT*T. 
S is then called a Kreiss Symmetrizer for the problem: 



(2.1) 



dxU = Au + /, {x > 0}, 



where / and g denotes respectively the Fourier-Laplace transforms of / and 
g; and u denotes the Fourier-Laplace transform of the solution u of the well- 
posed mixed hyperbolic problem (jl.ip . u is also solution, for all fixed C 7^ 
of the following equation: 



(2.2) 



Sdxii = SAu + SiAay^f, {x > 0}, 
Tu\x=o = ^9, 



Remark 2.1. Following our current assumptions, T is independent of C, 7^ 0, 
however, more general boundary conditions, of the form: 

r(C)u|x=o = r(c)5, 

can he treated. It would imply taking as boundary condition for (jl.ip : 
with for 7 big enough, 

where, T_{dt,dy) denotes the Fourier multiplier associated to that is to 

say is defined by: 

Referring for example to [3] and [7], Kreiss has proved that the exis- 
tence of a Kreiss symmetrizer for the symbolic equation is sufficient to prove 
the well-posedness of the associated pseudodifferential equation (here (jl.ip ). 
Indeed, multiplying by u and integrating by parts the equation: 

SdxU = SAu + S{Ad)-^f 

leads to the desired a priori estimates. For all (" 7^ 0, S'(C) is hermitian 
symmetric and definite positive on ker F. Let us sum up the properties crucial 
in the proof of the well-posedness of our problem: 

Proposition 2.2. For all ( = (r,7,??) such that + 7^ + Yl'jtiV'j = 1; 
there holds: 
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• S{Q is hermitian symmetric. 

• ^{SA) (C) := 1{SA + {SA)*){C) is positive definite. 

• —S{C) is definite negative on kerF and kerT is of same dimension as 
the number of negative eigenvalues in —S{C). 

Note that, by homogeneity of S, it is equivalent for the properties in Propo- 
sition [22] to hold for ICI = 1 or for |(^| > 0. As a consequence of the first 
point and third point of Proposition 12.21 the Lemma 11.41 applies and gives a 
matrix B such that: kerP = ¥,-^-{S^^ B). In the sequel, such a matrix B is 
fixed once for all. 

The following chapter contains a proof of Lemma 11.41 assorted of a detailed 
construction of B. 



2.2 Detailed proof of Lemma II. 4t Construction of the ma- 
trices B solving Lemma II. 4L 

As we will emphasize in next chapter, Lemma ll.4l is a crucial feature in our 
first method of Penalization. The aim of this chapter is to give a more com- 
plete proof rather than simply recalling Ranch's result and, in the process, 
to precise how the matrices B solving Lemma 11.41 are constructed. For all 
C 7^ 0, S'(C) is hermitian symmetric, uniformly definite positive on E_|_(^(^)), 
and uniformly definite negative on E^{A{()); as a consequence, S{C) keeps 
exactly p positive eigenvalues and N — p negative eigenvalues for all C 7^ 0- 
Basically, knowing that S is uniformly definite positive on kerP; we search 
to express ker P in a way involving S. Consider q € ker P, since, for all C 7^ 0, 
E_(5(C)) 0]E+(5(C)) = C^, we can split q in: 

q:=q+ + q' 
with q+ G E+(5(C)) and q- € E_(S(C)). 

Since dim/cerP = dimE+(S'((^)) = p, these two linear subspaces are in bijec- 
tion. Let us give the two main ideas behind this proof: one idea is to detail 
the bijection between q G kerT and q~^ G E+(5(C)) as it satisfies some con- 
straints, the other is to come down to the model case where the eigenvalues 
of S are either 1 or —1. Let us denote: 



-Mn-p 
Idp 



In a first step, we will prove the following result: 
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Proposition 2.3. If we assume thatY is a linear subspace ofC^ of dimen- 
sion p, and that there is C > such that, for all g G V, there holds: 

{S-\,q)>C{q,q), 

then the two following equivalent properties hold: 

• There is a hermitian symmetric, positive definite matrix R, such that: 

[qeY]^\!f^q(z E+{RSR) 
which is equivalent to: 

Y = E+{R^S). 

• There is a hermitian symmetric, positive matrix R, such that: 

[qeY]^ \Rq G E+{RS-^R)\ , 
which is equivalent to: 

Y = E+{S-^R^). 

Moreover, we can link the two properties by taking: 

R^ = SR^S. 

Proof. In this proof, we will show how to construct some matrices R 
satisfying the required properties. There is a {N — p) x p matrix of rank 
N — p such that < 1 and: 

Y={qe C^, q- = i<q+}, 

where g"'"[resp q~] denotes the projector on E+((5)~-'-) [resp E_((5)~^)]parallel 
to E_((5)-^) [resp E+((5)-i)]. Indeed, dimV = p = dimE+((^)-i), and 
= E_((5)-i)0E+((5)-i). Moreover, there ic C > such that, for all 
g G V, there holds: 



and thus 



{{S)-\,q) = -{q-,q-) + > C{q,q). 

- > C\q\^, 



11 



which imphes that < 1. We will show now that, for R constructed as 
follow: 



R 



Id„ 



there holds: 



Rq e E+{RS-^R) 



First, we see that the constructed R is trivially hcrmitian symmetric and 
positive definite since ||i^|| < 1. First, we have: 



RS-^R 



-IdN-p + NN' 






Idp - N*N 



and 



q- - ^q+ 
+ q+ 



Rq = 

Thus, since ||b^|| < 1, there holds: 

Rq € E+(^^-^^)j ^ [q- - Hq+ = O] ^ [g G V] . 
We will now prove that we have: 

{R)-^E+{RS-'^R) = E+{S-'^R^). 

Since RS~^R is hermitian symmetric, the linear subspace E^{RS~^R) is 
generated by the eigenvectors of RS~^R associated to positive eigenvalues. 

A basis of {R)-'^E+{RS-'^ R) is thus given by ((i?) ^Vj)j where Vj denotes 
an eigenvector of RS~^R associated to a positive eigenvalue \j. We have: 

RS~^Rvj = XjVj. 

Let us denote Wj = {R)~^Vj, we have then: 

RS"^ R?Wj = XjRwj ^ S~^R?Wj = XjWj. 

As a result, Wj is an eigenvector of S~^RF' associated to the eigenvalue Xj 
hence we obtain that: 

{R)-^E+{RS~'^R) = E+{S-^R^). 

We can also prove, the same way, that: 

RE+{RSR) = E+{R^S). 
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Now, taking 

Ix = SR S, 

we can check that: 

which concludes the proof. □ 
Lemma (??) is a Corollary of the following Proposition: 

Proposition 2.4. If denotes a smooth in(^ ^ 0, matrix-valued function 
in the space of hermitian symmetric matrices with p positive eigenvalues and 
N — p negative eigenvalues andkeiT denotes a linear subspace of dimension 
p and there is C > such that, for all q G kerF, there holds: 

{S-'q,q)>C{q,q), 

then the two following equivalent properties hold: 

• There is a smooth in (" ^ 0, matrix-valued function R, in the space of 
hermitian symmetric, positive matrices such that: 

[q G Kerr] ^ [VC / 0, R-\()q G E+{R{C)S{C)R{0)] , 
which is equivalent to: 

VC/0, Kerr = E+(^2(^)^(^))_ 



• There is a smooth in C ^ 0, matrix-valued function R, in the space of 
hermitian symmetric, positive matrices such that: 

[q G Kerr] ^ [VC / 0, R{Oq G K+{R{OS~HORiO)] , 
which is equivalent to: 

VC/0, KeTT = E+iS-\C)R\0)- 



Moreover, for all C 7^ 0, these two properties can be linked by taking: 

Proof. We will show here that Proposition 12.41 can be deduced from Propo- 
sition [2?3l Por all C 7^ 0) *S'(C) is a hermitian symmetric matrix, moreover S 
depends smoothly of C- As a consequence is also a hermitian symmetric 
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matrix depending smoothly of and as such, there is a nonsingular matrix 
V such that: 

5-1 = V* {s-^) V. 

Let us denote A the diagonahzed version of with eigenvalues sorted by 
increasing order, then there is Z depending smoothly of such that, for all 
C 7^ 0, we have: 

z*ic) = z~Ho, 

and 

A(c) = z*(c) i-s-') (cmc). 

As a consequence, V depends smoothly of ( since, for all C ^- 

where A is the diagonal matrix obtained by taking the absolute value of each 
eigenvalue of A. For the sake of simplicity, let us omit the dependence in 
Now, for all q G ker F, there is C > 0, such that: 

{S-'q,q) = {V*S~^Vq,q) = {S~\V q) , iV q)) > C {{V q) , {V q)) . 

Moreover V is nonsingular, thus there is C" > 0, such that, for all 
q ^ ker r, there holds: 

{S-'q,q)>C'{q,q). 

Moreover diml/^^ kerT = p, using Proposition 12.31 for all fixed ( ^ 0, there 
is a hermitian symmetric, positive definite matrix R{C)^ such that: 

y-i(C)kerr = E+((^(c))'^(C))=^(C)iE+(^(C)^(C)E)(C). 

We will now prove that we can construct R depending smoothly of C. First 
there is a (A'^ — p) x p matrix of rank N — p, depending smoothly of C, 
such that fore ah C / p(C)|| < 1 and: 

y-nC)kerr = {gGC^, = ^^(C)9+}, 

where g"'"[resp q~] denotes the projector on E_(_((S')~^) [resp E_((S')~^)]parallel 
to E_((S)-i) [resp E+((S)-i)]. R is given, for ah C / 0, by: 

^(C) = ^Js~Hc)R^{OS~HO, 
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with R given by: 

idN-p -m) ' 

Since S"^ = V* {S~^) V, there holds: S = V*SV, and, as a consequence: 

(VR)-^ ker r = E+{RV*SVR). 

As Ry*SVR is hcrmitian symmetric, a basis of the linear subspace K^{RV* SV R) 
is given by the eigenvectors of RV*SVR associated to positive eigenvalues. 
This leads us to consider vj = {VR)~^Uj satisfying: 

RV*SVRvj = XjVj. 

We have: 

RV*SVR{VR)-^Uj = Xj{VR)-^Uj. 

hence: 

{VR)RV*Suj = XjUj. 

Since {VR)RV* = {RV*)*{RV*) is hermitian symmetric and positive defi- 
nite, we can then define its square root. We define R by: 

R= \/iRV*)*{RV*). 

Since both R and V depends smoothly of (, so does R. Moreover, there 
holds: 

^Suj = XjUj, 

which gives: 

kerr = VRE+{RV*SVR) = E+{R^S). 

We have thus proved there is a smooth in ^ 7^ 0, matrix- valued function R, 
in the space of hermitian symmetric, positive matrices such that: 

[q e Kerr] ^ [VC 7^ 0, RT^Qq G E+(^(C)-S(C)^(C))] , 

which is equivalent to: 

VCt^O, Kerr = E+(^2^C)5(C)). 

Now consider R defined, for all C 7^ 0) by: 

R{0 = Vs{C){R{0?s{0, 



m = 
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^(C) = \/{R{C)v*{C)Sic)ym)v*{c)S{C)). 

C 1-^ -R(C) is smooth and, for all R{C) is a hermitian symmetric, positive 
definite matrix. Moreover, there holds: 

[q G Kerr] ^ [VC ^ 0, R{Oq G E+(i2(C)5-i(C)i?(C))] , 

which is equivalent to: 

VCt^O, Kerr = E+(S-i(C)i?'(C)). 

Let us detail the computation of -R(C)- 

R{0 = \/s{Ov{Or\ov*{OS{0. 

Moreover 

r'{0 = s-\OR'{Os-\o, 

we have thus: 

^(C) = ^(i2(C)^-HC)v^*(C)5(C))*(i?(C)5-nC)^*(C)5(C)), 

which gives: 

5(C) = (i2(C)5-nC)^*(C)^(C))* (i?(C)^-HC)v^*(C)^(C)) • 

We recall that R is given, for all C 7^ 0, by: 
and that for all C 7^ 0; ^(C) is given by: 

v{0 = (A(C))-^^(C), 

where 

A(C) = z*{0 i-S'') (C)z(c) 

with A is a diagonal matrix with real coefficients: (Ai, . . . , Xn), and A de- 
notes the diagonal matrix with diagonal coefficients (|Ai|, . . . , |Ajv|). 

Remark 2.5. In the construction of B the only freedom we have resides in 
the choice of 'i^. 

□ 
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2.3 A change of dependent variables. 

Let us denote by R := B2 and v := Ru. v is hence solution of ()2.3p : 

(2.3) 



R-^SR'^d^v = R-^SAR'H + R-^S{Ady^f, {x > 0}, 



[ TR ^v\^=o = Tg, 

We will adopt the following notations: Sr := R^^SR^^, Ar := RAR^^, 
and Tfj := TR^^. We first observe that: 

kerVn = RkeiT = RK+dSy^ R^). 

but Sj^^ = RS-^R thus 

kevTR = RE+{R-^SbR) = E+{Sr). 

This is where Lemma 11.41 is used in a crucial manner. Let us denote by 
the projector on E_(S'ij) parallel to K^{Sr) and by by the projector 
on E+(S'/j) parallel to E_(S'ij); P~ and P"*" denoting the associated Fourier 
multiplier. Since Sr is hermitian symmetric, P~ is in fact the orthogonal 
projector on E_{Sr). The problem (|2.3|) can then be written: 

f SrO^v = SrArv + R-^S{Ad)-^f, {x > 0}, 

\ p-i)U=o = p-r^, 

This problem is well-posed because, as a direct Corollary of Proposition 12.21 
we have: 

Proposition 2.6. For all ( such that + 7^ + \r]\'^ = 1, there holds: 

• Sr{Q is hermitian symmetric. 

• ^{SrAr) (C) is positive definite. 

• —Sr{C) is definite negative on kerL/j and the dimension ofkeiTR is 
the same as the number of negative eigenvalues of —Sr{C). 

Proof. For the sake of simplicity, let us omit the dependence in in our 
notations. 

• Sr := R~^ S R~^ , and both S and R are hermitian thus Sr is hermi- 
tian. 
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• SrAr = R-^SAR-^, thus for all q G C^, there holds: 

2mSRAR)q,q) = {SRARq,q) + {q,SRARq) = {R'^ SAR-\, q) + {q, R'^ SAR' 

since R~^ is hermitian, we have then: 

= {SAR-\,R^\) + {R-\,SAR-\) = 2{^{SA)R-\,R-\). 

Since ^{SA) is positive definite and R is invertible, 3? (SrAr) is thus 
positive definite. 

• By construction of R, it satisfies kerF^ = E_|_(S'ij,), with Sr hermitian. 
As a consequence —Sr is definite negative on ker Tr and the dimension 
of kerF/j is the same as the number of negative eigenvalues of —Sr. 

□ Let us mention that, since R and S remains uniformly bounded in 
C / 0, / and R~^S{Aci)~^f belongs to the same space. In a same spirit as 
[5], this suggests the following singular perturbation of (|2.3p : 

SRd^v' - ^F-v'l^^o = SrArv' - ^p-r5l,.<o + R-^S{Adr^f, {x e M}, 
This is equivalent to perturb ()2.2p as follow: 

Sd^u' - ^RF-Run,^o = SAu' - ^RF-Tgl,^o + {x € M}, 

Finally, this induces the following perturbation for (jl.lj) : 




nu' + -Mn=l,<o = /l,.>o + -^l..<o, {x G M}, 

e e 

u^\t<o = 0, 



where 

M := -e^^AdS-^RF-Re-''\ 

e = -e^^AdSr^RF-T^, 

and S_{dt,dy) [resp denotes the Fourier multiplier associated to 

5(C) [resp RiC)]. 

3 Proof of Theorem 11.61 

First, we construct an approximate solution of equation (|2.4p (which is also 
equation (jl.4p ) , then prove suitable energy estimates that ensures yf and its 
approximate solution both converges towards the same limit as 
e^0+. 
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3.1 Construction of the approximate solution. 

is the solution of the well-posed Cauchy problem: 

Hyt + ^Mn^l^<o = /la;>o + ^^la:<o, {x G M}, 

is moreover the solution of the well-posed Cauchy problem: 

^ SA-^Uyt + -^SA-^Mun:c<Q = SA'^f l^>o + ^SA-^ei^<o, {x e M}, 
^lf\t<o = 0. 

The associated equation after tangential Fourier-Laplace transform writes : 
Sd^u' - -BF-Ru'l^<o - SAu' = --RF-Vgl^^o + S{Aa)-^ fla,>o, {x G 

£ £ 

or alternatively: 

SrO^v' + -F-v'l^^o = SrArv' + -P-r5l^<o + R-^S{Aa)-^f, {xeR 

£ £ 

We will use the following formulation as a transmission problem in our con- 
struction of an approximate solution: 

SRd^v'+ = SrArv'+ + R-^S{Ad)-^f, {x > 0}, 

SrO^v'- + -F'v'- = SrArv'- + -p-r^, {x < 0}, 

£ £ 

For n an open regular subset of W^'^^, and p G N, let us introduce the 
weighted spaces H^{Cl) defined by: 

H^{n) = {zue e^'L\^), ||ti7||^.(n) < oo}; 

where 



!2 

a,\a\<e 



We will construct an approximate solution v^pp of 'if. y^pp will be con- 
structed as follow: 
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where is an approximate solution of n^^ satisfying the foUowing ansatz: 



M 

e± 



ILapp 



j=0 

I — 7 _|_ _|_ J 

where the profiles C/j" belong to i^^ ^ (O^), where stands for [0, T] xR^. 
Denote 

is then an approximate solution of v''^ and is of the form: 

M 



'^app 



j=0 



where 

and conversely 



The profiles U_j can be constructed inductively at any order. Let us show 



-J 

how the first profiles are constructed: Identifying the terms in gives: 

Hence, P^V^^ remains to be computed in order to obtain the profile 

Identifying the terms in gives then that V^^ is solution of the well-posed 
problem: 



(3.1) 



SRd,V+ = SnARV+ + R-^SiA^r^f, {x > 0}, 



The associated profile 



belongs then to H^i^^)- ^oi' 

eover, the problem (j3.ip is Kreiss-Symmetrizable 
and thus the trace of the profile U_q , see [3j for instance, satisfies: 



I/+U=oG^'(Tt) 
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Since V^^ has just be computed, Vq |a;=o is given by: lo^|a;=o — U=o = 
and thus, there holds: 

P V(^\x=0=^~Vq \x=0- 

Moreover 

SrOxV^- - P-Ff = SrArV^-, {x < 0}. 
Projecting this equation on E+(S'r) coUinearly to E.^{Sfi) gives then: 

SrOxF+V^- - F+SrArV^- = 0, {x < 0}, 

Since 

F+SrArV^- = F+SrArF+V^- +F+SRAR¥-rg, 
we have then: 

SRdxiF+V^-) - F+SrAr{F+Vq-) = F+SRARF-Tg, {x < 0}, 
and as a consequence, P^V^^ is solution of the following problem: 

(SRdx{F+VQ')-F+SRAR{F+Vo-)=F+SRARF-Tg {x<0}, 
\P+yo"U=o = P+VU'=o- 

Let us precise how (|3.2p has to be interpreted: we denote w = F^V^^ . w is 
then totally polarized on E_|_(S'/j), and satisfies the problem: 

F^w = w 

(3.3) I SrOxW - F+SrArw = F+SRARF-Tg {x < 0}, 

^u;U=o = P+^o^U=o. 

As we will see, the problem (|3.3|) is Kreiss-Symmetrizable and thus well- 
posed. Indeed, for all C such that + 7^ + |r/p = 1, we have, omitting the 
dependencies in in our notations: 

• For all q G C^, there holds: 

{SRq,q) = {q,SRq). 

• Since Re{SRAR) is positive definite and P+ is an orthogonal projector, 
there is C > such that, for all q G E+(S'j?), there holds: 

{F+SRARF+q,q) + {q,F+ SRARF+q) > C{q,q). 

Indeed, for all q G E_|_(S'/j), there holds: 

{F+SRARF+q,q) = (P+^^^flP+g, P+g) = {SRARF+q,F+q). 
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• —Sr is definite negative on kerP+ that is to say, that there is c > 
such that, for all q G kerP+, there holds: 

{-SRq,q) < -c{q,q). 

Moreover kerP"*" has the same dimension as the number of negative 
eigenvalues in —Sr. 

The profile U_q can then be computed by: 

[/- := e^*j?^-^ {R-^{w + F-Tg)) 

belongs to H!^{CI^), moreover its trace U^\x=o belongs 
to H!^{Tt). Consider now the equation: 

^-Vf = SrO^Vo- - SrArV^- , {x < 0}. 
Since F~V^\x=o = 'F~V^\x=o, is solution of the well-posed problem: 

f SRdxV+ = SrArV+, {x > 0}, 

\ P-V+U=o = SRdxVo-\x=o - SrArVo-U=o- 
Due to the loss of regularity in the boundary condition, the associated profile 

belongs to ^(fij), moreover its trace U_f\x=o belongs 
to Hj ^(Tr). Moreover, applying P+ to the equation: 

F'V^- + SrArV,- = SRdxV{- , {x < 0}, 

we obtain: 

f SBd,F+Vf = F+SrArF+V{- + F+SrArF-V^, {x < 0}, 
\f+V^-\x=o = F+V+\x=o. 

As before, let us take P+Fj^~ as the unknown of the well-posed problem: 

f SRdx{F+V,-) - F+SrAr{F+V^-) = F+SrAr {SrOxYq- - SrArV^-) , {x < 
\{F+vn\x=o = P^V+\x=o. 

This problem is Kreiss-Symmetrizable since, for all ^ such that 
+ 7^ + |?7p = 1, there holds: 
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• For all qeC^, there holds: 

{SRq,q) = {q,SRq). 

• There is C > such that for all q G E-|-(5'i?), there holds: 

{F+SRARF+q,q) + {q,¥+ SRARF+q) > C{q,q). 

• —Sr is definite negative on kerP"*" that is to say, that there is c > 
such that, for all q G kerP+, there holds: 

{-SRq,q) < -c{q,q). 

Moreover kerP^ has the same dimension as the number of negative 

eigenvalues in —Sr. 

However, due to a loss of regularity in both the source term and the boundary 
condition, the associated profile 

U- = e^T-^ (P+l/f + SrB^V^ - SrArV^)) 

belongs to i/^ The construction of the following profiles can be 

pursued at any order the same way. In practice, we take: 

As a result, the approximate solution writes u^^p := v^pplx>o + Uapp^x<o', 

where v^pp belongs to Hi^{0,^) and v^pp belongs to ^(il^). 
solution of a well-posed problem of the form: 



(3.4) 



Hv^pp + -Mi^ppl^^o = /l^>o + -eix<o + erf, {x G M}, 



I, Mapp |t<0 = • 



Where := r^+lx>o + z:^"lx<o, with € iJ^' ^ and 

Remark 3.1. In the case where g = 0, the loss of regularity in the profiles 
is delayed by one step. As a result, in this case we obtain: 

y^pp e H!;{n-), 
r'+ e H!;in+), 



3 
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3.2 Stability estimates 

We will begin by proving energy estimates on the following equation: 
(3.5) SRAne' - SrO^c' + ^P-e^l.<o = ef', {x G M}, 



where = i? {J^{e ^^u^) — T{e ^^Uapp)) ■— with w'^ = — 



Refering to ()3.4p . is the solution of the Cauchy problem: 



+ -Mu;^l^<o = erf, 
(3.6) { e 

For a fixed positive e, the perturbation is nonsingular and thus the principal 
part of the pseudodifferential operator 7i + is the same as the principal 
part of Ti.. Hence, there is a unique solution of the Cauchy problem (13. 6p : 
which belongs to H^~^{Q.t)- In order to simplify the notations, in this 
chapter we shall denote by and HS/ the spaces: L?'{Q.t) and H^{Q.t)- 
We recall the definition of the weighted spaces: H^{VtT) for p € N. 

where 

ct,\ce\<g 

For fixed positive e, there holds: 

/oo 
-oo 



-oo 

fOO 



/oo 
-oo 

Using the equation, we have then: 

r°° 1 

/ Re{SRARe^ + -P~e^ - er^ , e^)i2 dx = 0. 

J —oo ^ 

which is equivalent to: 

/oo -j^ /"OO 
Re{SRARe^,e')L2 dx = — / i?e(P~e^ef% e^)i2 dx 
-oo ^ J —oo 

/oo 
Re{r'^ , (f) dx. 
-oo 
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(<^-|;)7rfe+Jiip-riii.<i^irr.. 



But Re{SRARe^,eF) = {Re{SRAn) eF ,e^) and Re{SRAR) is positive defi- 
nite, hence there is C > 0, independent of e such that: 

Cl\W\\l2 + - j i^e(p-e^e") < / Re{er' , e^) ^2 dx. 

Thus, because P~ is an orthogonal projector, for all positive A, there holds: 

C^Trili. + ^||P-e^|li. < ^^lle^li. + ^llsf^lli. 

Choosing A big enough we have C—j^ > and the following energy estimate: 

-.2- 

This shows that converges towards zero in when e tends to zero, with 
a rate in 0(e). We recall that is given by: 

:=RJ^{e-'y\y^pp-u')), 

and is given by: 

:= i?JP (e-T*r^) . 

Moreover, since R and P~ are two uniformly bounded, uniformly definite 
positive matrices, there are two positive real numbers a and P such that, 
for all C 7^ and x eM., there holds: 

. a\\F-T{e-''\ul^^-u^)) < \\F-e^\\l,. 



\r 



e\\2 



Applying then Plancherel's equality we obtain then: 

(c - ^) ilKpp - u'\\l-rtL2 + ^||M {y^^ - 2f) \\l,tL2 < ^'-^Wr.ll.tL^- 

We have thus proved there are two positive constants c and C such that: 

1 Cs"^ 

ClWulpp-u'Wl-ytL^ + -||M {V^pp-U') \\l,tL2 < —\\f\\l,tL2. 



Let us denote by ||.||^, := + Un^^^-y More generally, when 

G H^, there is two positive constants Cp and Cp such that: 



— 1 1 jl^ 



As we have seen during the construction of the profiles, ^ = A; — 3 in general 

3 
2 



and p = A; — I in the case where g = 
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3.3 End of the proof of Theorem 11.61 

As a consequence of our stability estimate, there holds: 

Moreover, by construction of u^pp, there holds: 

Hence, we have proved that: 

Wn^ - u llHfc-3(c-) + 11-^ ~ '"IIh'=-3(q+) = C'(e^). 
By the same arguments, if g = 0, there holds: 

\\u'-u-\\\ 3 +\\u'-u\\\ 3 ^ =C(e')- 

This completes the proof of Theorem ll.6[ 

4 Proof of Theorem 11.91 

Like in the proof of Theorem 11.61 we begin by constructing formally an 
approximate solution of equation (jl.7p . We prove then suitable energy esti- 
mates that ensures both and its approximate solution converges towards 
tt as e ^ 0"'". 

4.1 Construction of an approximate solution. 

The goal of this Lemma is to replace the boundary condition rn|2:=o = of 
problem by a condition of the form 'P_~ i^~'''^u)\x=o = h with a suitable 
h G H''{Tt). 

Lemma 4.1. Let u denote the unique solution in H^{Q,^) of the mixed 
hyperbolic problem (ll.ip . F.~^{dt,dy,j) [e'^'^^uj does not depend of the choice 
of the boundary operator T and of g. Let us introduce the function h of 
H''{Tt) defined by: 

P- (e-^*^;|,=o) + n {e-^\g - v\,=o)) ■ 
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The solution u of the mixed hyperbolic problem (jl.ip is the unique solution 
of the following well-posed mixed hyperbolic problem (|4.ip ; 



(4.1) 



'nu = f, {x> 0}, 

, u\t<o = . 



In addition, the mapping (/, g) —* h is linear continuous from 
H^{Q+) X H^{Tt) to H^{Tt). 

Proof. Let v denote a solution in H'^{Q.t) of the equation: 

( Hv = f, {t,y,x) e ^T, 
\v\t<o = . 

We introduce then U which is, by definition, the solution of the following 
mixed hyperbolic problem: 

nv = 0, {x> 0}, 

T{dt, 7)UU=o = Hidt, dy,-f)g - T{dt, dy,-f)v\x=o, 
U|t<o = . 

The right hand side of the boundary condition is, a priori, 

in H^~^{Tt)- Hence the solution U belongs to H^^2{n)p). By construction 

we have: 

(4.2) u = V + v. 

Hence, since u € H^{^^) and v G H''{^^), in fact we have: 

Let U denote the Fourier-Laplace transform in (t, y) of U (Fourier-Laplace 
transform tangential to the boundary) given by: J^(e-^*U). It satisfies the 
following symbolic equation: 

r d^i] = ^(c)u, {x > 0}, 

\ r(C)uU=o = nOg - r(C)^U=o, 

where g and v denotes respectively the tangential Fourier-Laplace transform 
of g and v. Since A{C) is independent of x, projecting the above equation on 
E+(A(C)) gives then: 

a^p+U = A{C)P+t. 
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Moreover P+lLJ|3;=o G E_{A{C)) f}E+{A{()) since lim^^oo P+U = 0. Hence, 
there holds: 

P+U = 0, 

and thus 

u = p-u. 

The boundary condition: 

r(C)U|^=o = r(c)5-r(C)«|x=o 

is equivalent to: 

U|x=o e 5 - v\x=o + Kerr. 

We have thus: 

P~UU=o e 5 - v\x=o + ker T. 

Let us denote by 11 the projector on E_(^) parallel to kerT, which has a 

sense because the Uniform Lopatinski Condition holds. 

Since Ul2:=o € E_(A), and of the Uniform Lopatinski Condition, the above 

boundary condition is equivalent to: 

np-u|^=o = n(5-i)U=o), 

and thus, because P~U|j;=o belongs to E_(^), we have: 

p-lLJU=o = n(5-{;|^=o). 
As a consequence, we obtain: 

P~u\x=o = P~*U=o + n(^ - v\x=o)- 

Hence, there holds: 

P- (e-T*uU=o) = P- {e-^'v\x=o) +n{e-^\g - v\x=o)) := h. 

P_'^{dt, dy,j) [e~'^^u) = ^{dt,dy, 7) (e^^*u) , thus it does not depend of the 
choice of the boundary operator V and of g. Moreover, since u\x=o G H'^{Tt), 
it follows that 

g G H^{Tt)- Now, since the Uniform Lopatinski Condition holds, u satisfies 
the following energy estimate: 
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More generally, we have: 



1 

7' 



u 



h = P~(e~'^*n|2^=o) hence 



where \\w\\l, := EL=o t'"!"' Il<9"inii2 



L2{Tt) 



u\x=o\\\2(j^) — C||'wU=o|le7tL2(Xy); 



and for < j < d — 1, there holds: 

||5,/i|li.(x,)<c,-||r?,P-^(e-^^ 
More generally, we have: 



=o|| < Cj-|ln|2,=o||Hi 



But 7 is a positive real number fixed once and for all at the beginning of the 
paper, hence this proves that the mapping (/, g) ^ h \s continuous from 

X H^{Tt) to H^{Tt). □ 
As we will see. Lemma l4.1l is central in our construction of an approximate 
solution. We will construct an approximate solution 



"app 



along the following ansatz: 



M 



"app •~ / > "-^^7 {^jyj^)i 



j=0 



withn+G/7^ 2^(0+), M+|^.=o G -??7 ' (Tt); and 



"app 



M 
j=0 



Uj {t,y,x), 



with u~ G if-y ^ (^2^), u^\x=o £ ^ {'^t)- As usual, we will refer to the 
terms u- as profiles. We will rather work on the reformulation of 
problem p.7p as the transmission problem (j4.3p : 

r7^^x^+ = /, {x>0}. 



(4.3) 



1 



^rfe^*/i, < 0}, 



x=0 - 



=0 = 0, 
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Plugging u^+p and u^app (|4.3p and identifying the terms with same power 
in e, we obtain the following profiles equations: 

• Identification of the terms of order : 

(4.4) ^rfeT*p-e-^*no = Ade^% {x < 0}. 

• Identification of the terms of order e*^ : 

(4.5) Huq + Ade'*¥re~^\^ =0, {x < 0}. 

(4.6) Hut = /, {x > 0}, 

• Identification of the terms of order for j > 1 : 

(4.7) Uu] + Ade^*p-e~^*n7^i = 0, {x < 0}. 

(4.8) Uu^ = 0, {x> 0}, 

• Translation of the continuity condition over the boundary on the pro- 
files: 

For all 1 < j < M, there holds: 

(4.9) it+|^=o -'"7U=o = 0. 

Denote by := T{e~'^^u^) . We have then: 

We will now give the equations satisfied by the Fourier-Laplace transform 
of the profiles: . First, equation (j4.4p is equivalent to: 

p-iig = T{h)^ {x < 0}. 
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We deduce from this equation that there holds: 

Then, using (|4.9p for j = 0, and (|4.6|) gives that, for 7 big enough, 

<=^(e-^*n+), 
where Uq is the solution of the well-posed first order ODE in x: 
(d,u^-Au+ = T{e-''\Ad)-'f), {x>0}, 

Thus Uq is solution of: 

r nu+ = /, {x > 0}, 

\e^*p-e-^*<U=o = /i. 

Thanks to Lemma [4. 11 we recognize Uq as the solution of our starting mixed 
hyperbolic problem p.l|) . Once Uq is known, so is Uq and thus u^\^^q is 
given by: 

Uq \x=o — ^0 U=o- 

Moreover, 

Uq\x=o = U'=o £ H^i^r)- 

By there holds: 



dxUQ — Auq + P -U]^ = 0, {x < 0}. 
As a consequence, F^u^ is given by the well-posed ODE: 

f 5,(P+^io ) - ^(P+^i^) = 0, {x < 0}, 

I P + Ug j^ = = P^lio U' = 0- 

Indeed, since kerP+(C) = E_(j4(C)), this problem satisfies the Uniform 
Lopatinski Condition: for all C 7^ 0; there holds: 



E_(A(C))0E+(^(C))=C 



N 



For 7 big enough, by linearity of the inverse Fourier transform, Uq can then 
be computed by: 



- := e^*^-i(p-no ) + e^*^-i(P+no ). 
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Following up with that process of construction, we can go on with the con- 
struction of the profiles at any order. Indeed, assume that all the profiles 
{u~j,uj) up to order j have been computed. Then consider the equation 
obtained through identification: 

^"""7+1 = -dxuj + AuJ, {x < 0}. 

We see there is a loss of regularity between u~_^-^ and uj . 

Let us say that G Considering the traces, we have then: 

uf\x=o G ^(Tt)- We will show in this part how the Sobolev regularity 

of the profiles Uj^-^, which is by definition rrij+i, can be computed know- 
ing rrij. To begin with P_''uJ_^_^ belongs to H"''^ ^(0,^). 'PrUj'^-^\x=o, which 

belongs to iJ™' ^{'^t), is known by P~uj^^\x=o = e^* J^'^ (P' u'^^^\x=o) , 
with: 

P~«/+iIj:=0 = P"'"7+i|x=0- 

Hence, Uj^j^i '■= ^(e~'''*n^^) is the solution of the first order ODE in x : 

f t^x-u++i - = 0, {x> 0}, 

\ P~^i^+iU=0 = 'P'uj^i\x=0- 

Since kerP~(^) = ¥,^{A{0), this problem satisfies the Uniform Lopatinski 
Condition: for all ^ / 0, there holds: 

E_(^(C))0E+(^(C)) = C^. 

3 

As a consequence, this problem is well-posed and, uj^^eH^' 2 More- 
over, there holds: 

3 

u^+i\x=o = uj+i \x=o e -H"^' ^ (Tt)- 

Indeed, P+ii++i G hence P+u++J^=o £ //^'"^(Tt) and thus 

u'j_^_i\x=o € -ff™^ ^C^t)- Furthermore, we have: 

"f^j-l-l U=0 ~ '^j^4-iU=0- 
Applying P+ on the following equation: 

P~«7+2 = -dxuj+i + AuJ_^-^, {x < 0}; 
we obtain then the equation: 

dx{P+uJ+i) - AP+uJ^^ = 0, {x < 0}. 
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Remark 4.2. 

shows that the "Fourier profile" fij+i must be so that —dxuj^i + -^uj+i is 
polarized on E_(j4). It is indeed the case because we search for uj^i satisfy- 
ing: 

dx{P+uJ+,)-AP+u-^^ = 0, {x<0}. 



u 



'j+i is given by: 



u 



3 

with P+wJ+i = e^*^-i(P+n-^J belongs to '~^{^t) and is the unique 
solution of the well-posed first order ODE: 

r 5,(p+u-+i) - ^(p+uj+i) = 0, {x< 0}, 

The profile uj^-^ belongs to Hj ^ '^(Q^). This achieves to show that the 
knowledge of {uj',uj), allows us to compute {uj'_^_i, uj_^_^) . 
Moreover m^+i = rrij — |, that is to say that a construction of each supple- 
mentary profile consummate | of Sobolev regularity. In practice, we take: 

,,£+ _ ,,+ 
"'app ~ "O ' 

As a result, the approximate solution writes u^pp := u^^plx>o + u^pplx<o', 

where belongs to H^{^1^) and u^pp belongs to ^ (O^). The so de- 
fined u^pp is solution of a well-posed problem of the form: 
(4.10) 

Hulpp + ^Ade^'p-e-"''ulpplx<o = flx>Q + ^Ade^'hlx<Q + £r^ 

■"app |t<0 = . 

Where r^ := r^^la;>o + r^^^xKO, with € ^ (^t) 

r'- e Hl;-^{n:^). 
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4.2 Asymptotic Stability of the problem as e tends towards 
zero. 



Denote by 



u 



app 



. By construction of u^^p, v'^ is solution of the 



following Cauchy problem: 



1 



(4.11) 



e 



For all positive e, this problem is well-posed. In order to investigate the 
stability of this problem as e goes to zero, we will reformulate it as a trans- 
mission problem. The restrictions of to {x > 0} and {x < 0}, respectively 
denoted by v^^ and v'^~ are solution the following transmission problem: 



(4.12) 



' nv^+ = er^+, {x > 0}, 
nv'~ + ^Ade^^p-e~^^v'- = er^~, {x < 0}, 



U^"^|x=0-f^ U=0 



0, 



Let us denote by the function, valued in M."^^ , defined for all {x > 0} 
and (t, y) G [0, T] x M'^-i by: 



V'it,y,x) 



V'+{t,y,x) 
V'-it,y,-x) 



is solution of the Cauchy problem (j4.1ip iff is solution of the mixed 
hyperbolic problem on a half space (j4.13p given below: 



(4.13) 



where 



f + B^V = eR\ {x > 0}, 

t<0 



-1 



n = dt + Y, 



A, )^^ + 











Mde^*P 



Ad 


--yt 





-Ad 



dx, 



R'{t,y,x) 



r^+{t,y,x) 
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and 

r={ld -Id ) . 
Returning to the construction of our approximate solution, we have 
R'' G H^~^n+) X H^-''^{n+) and is such that R'\t<o = 0. 
In fact i?^ G HI^' (Q^) with k' = k — 3. For all positive e, there exists a unique 
solution in H^^^{fl^) to the above problem. We will prove here that this 
solution converges, uniformly in e, towards in (^^-r), as e vanishes. As 
in the proof of Kreiss Theorem, see p] for instance, existence of solution 
for mixed hyperbolic systems like (jl.7p or (j4.13p . are obtained through the 
proof of both direct and "dual" a priori estimates on an adjoint problem. This 
estimates results in the constant coefficient case of estimates on the Fourier- 
Laplace transform of the solution. Additionally, if this "Fourier" estimate 
can be proved, both direct and "dual" energy estimates are deduced from 
it. In a first step, let us recall formally how to conduct the Fourier-Laplace 
transform of a mixed hyperbolic problem: 

'nu = f, {x> 0}, 

< Tu\cc=o = 9, 
_ u\t<o = , 

Denote by := e~'*'*n, is in particular a solution of the following problem: 

Hu^+-iu^ = e-^^f, {x>0}, 

ru*|a;=0=5' • 

We take then the tangential (with respect to (t,y)) Fourier transform of this 
equation, which gives: 

d-i 

AddxU* + (7 + «t)u* + ir]j ^ AjU^ = T {e~^'^ f) , {x > 0}, 

i=i 

^r-u*|j.=o = 5 • 

Multiplying this equation by A^^ , we obtain that u* is solution of the fol- 
lowing ODE in x: 

f d^u, - Au, = {Ad)-^J^ (e-^V) , {x > 0}, 

I ru*|a;=o = 9 ■ 
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Note that, u^, and u can be freely deduced from each other through the 
formulas: 

and 

We shall now introduce a rescaled solution Vf of the solution of 
()4.13p defined as follows: Vf{t, y, x) := V^{t, y, ex), and the rescaled remain- 
der: R^{t,y,x) := R'^{t,y,ex). Denoting hy y_ = J^{e'^^y_), the associated 
equation writes then: 



dxV - eAV + MV = e^i^^ {x > 0}, 



where 



We remark that 



= I S P "(C) 



5.4(C) = AH) = A{0. 

with C = (f,7,57) := e^. Moreover P~ is homogeneous of order zero in ^. 
Let us denote R^{C,x) := R^{(,x) and F (C)^:^) •= Y. (Ci^;). Hence F is 
solution of the following problem: 



-A(C)+M(C) V =e'R'iC,x), {x>0}, 



fZ%=o = . 

As a consequence, the Uniform Lopatinski Condition for problem (|4.13p 
writes: For all 7 > 0, 

|(iet(E_(i(C) -M((:),kerr)[ > C7 > 0. 

In view of the proof of the Proposition (j4.3p . we recall that the spaces E±(A) 
have to be considered in the extended sense defined above. 

Proposition 4.3. Since Ti satisfies the hyperbolicity Assumption in As- 
sumption \l.li the Uniform Lopatinski Condition is satisfied for our present 
problem; that is to say that, for all C such that 7 > there holds: 

|det(E_(i((f) - M{C),keTr)\ > C >0. 
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Proof. We will begin to show that the Uniform Lopatinski Condition writes 
as well that for all C ^ there holds: 

(4.14) E+(A(C)-P-(C))nJE-(^(C)) = {0} . 

This notation keeps a sense for such that 7 = because we will prove 
a posteriori that the involved linear subspaces continuously extends from 
{C,7 > 0} to {C)7 = 0}. Then we will prove that, for all C, the property 
14.141 holds true. The Uniform Lopatinski Condition writes actually, for all 

E_(i(C)-M(C))nkerf = {0}. 
and thus, since we have: 

E_(i(C) - M(C)) = E_(^(C)) X E+(A(C) - P-(C)), 

and by definition of F, the Uniform Lopatinski Condition writes then that, 
for all (" 7^ 0, there holds: 

E+(^(c)-p-(c))n^E-(^(c))={o}. 

Lemma 4.4. 

E_ fA(C)-P-(C)) =E_ (a{0 



E+ (A(C)-P-(C)j =E+ (^(C)j. 

Proof. For all (" 7^ 0, there is an invertible N x N matrix with complex 
coefficients P{C) such that: P~^AP is trigonal and the diagonal coefficients 
are sorted by increasing order of their real parts. Let us denote by v the 
dimension of E_ (A) . The above matrix P traduces the change of basis from 
the canonical basis of into (f 1, . . . , Ujy, v,yj^i, . . . , V]y), where 

Span {ivk)i<k<u) = IE- (A) , 

and 

Span ((ffc)i.+i<fc<Ar) = (A) . 

Moreover, there holds 

p-ip-p = D 

where D is the diagonal matrix whose first diagonal terms are equal to 1 
and the N — v last diagonal terms are null. 

P-^{A - P-)P = P-^AP - D. 
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P AP—D is also trigonal, with the same eigenvalues with positive real part 
as P~^AP and the same eigenvalues with negative real part as P^^AP — Id. 
As a consequence, for all C 7^ 0) there holds: 

E_ (a(C)-P-(C)) =E_ (^(O), 
E+ (a(C)-P-(C)) =E+ (^(O). 

□ 

As a consequence of Lemma 14.41 the rescaled Uniform Lopatinski Con- 
dition for e > 0, e — > happens to be exactly the same as the one written 
for bigger positive e. Indeed, it writes, for all 7^ : 

E+(^(C))nE-(^(C)) = {0}. 

□ The Lopatinski condition is satisfied, and, as a result, the following, 
uniform in e, energy estimate holds for Vf , for all 7 > 7a;' > : 

^,||T/e||2 I ||T/e|2 11 ^^llc-P£||2 

which is equivalent to: 

(4-15) 7ll^^ll^.'(o.) + II^^U=o||^..(^^) < ^\\eR%,^^,^ . 

This proves the convergence of towards zero in H'^'{n+). The weight 
7 is fixed beforehand thus, in fact, the solution of ()4.13l) tends to zero in 
(0.^) at a rate at least in 0{e). 



5 End of proof of Theorem 11.91 

Let us consider defined by: 

This notation is perfectly fine because the so-defined function is solution of 
an equation of the form (j4.13p . Moreover, thanks to the stability estimate 
(|4.15p . there is 7^ positive such that, for all 7 > 7^,, we have: 

2 2 2^2 

7ll<pp-^^^ll^fc-3(f^+)+7ll<pp-'"^ll^fe-3(f^-) + ||<pp-^i1l^fc-3(^^ < — Ik-R^ 
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Hence, it follows that: 

,e||2 



,e||2 



Moreover, by construction of n^^^, we have: 



"app 



app 



U 



As a result, we obtain that there holds: 



u 



Oie' 



This concludes the proof of Theorem II. 9[ 



6 Appendix: answer to a question asked in 

In this chapter, we will show that the loss of convergence observed numeri- 
cally in [11] in a neighborhood of the boundary is due to a boundary layer 
phenomenon. We consider the 1-D wave equation: 

{ duU -c^dxxU = G]0,7r[xM+, 

f/jx = = U\x=TV = 0, 

U\t=o{x) = sin{x), 
I dtU\t=o = 0. 



(6.1) 



As in [TT], we define then = C/^+l^>o + U^~lx<o by: 

) G]0,7 

0, g] - oo,0[xl 



(6.2) 



duU'+ - c^dxxU'^ = 0, (x, t) G]0, ^[xl 



dttU' 



1 



^^dxxU'- + -^U'- 



: 



\x=0 







C/^ |a;=0 — \x=0 

dxU'+\x=o-dxU'- 
U'-^\x=n = 0. 

t=o[x) = sin{x), {±x > 0}. 



Jje±\ 



\t=o = 0, {±x > 0}. 
in 

fying the following ansatz: 



We will now construct formally an approximate solution U^^p of U^^ satis- 



^ app 



M 

Y,U;{t,x)e\ 
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M 



^ app 



j=0 



where the profiles Uj (t, x, z) := U_j (t, x) + U* (t, z), with 



lim e^^'^C/r 



0, 



for some a > 0. Since the stability estimates are trivial here, we will only 
focus on the construction of 

^app ^app^xyO + Ui^pplx<0- 

Plugging U^^p into problem ()6.2p and identifying the terms with same power 
of e, we obtain the following equations: 



Moreover, Uq~ = as it is the only solution of the problem: 



u^- - c'd,,u*- = 0, {z<0}, 

lim U*- = 0. 



^app converges towards Uq as e ^ 0+. As awaited Uq is the solution of the 
well-posed 1-D wave equation: 

' dttU^ - c'dx.U^ = 0, (x, t) G]0, 7r[xM+, 
Uq\x=0 = U_o \x=0 + \z=0 = 0. 



0. 



\X = TT 

U^\t=o{x) = sin{x), {x > 0}. 
[dtU+\t=o = 0, {x>0}. 

Let us write the following profiles equations: First, we can see that, for all 
j > 1, there holds: 



0. 



where is the solution of the well-posed profile equation: 



-dttU*- = 0, {z< 0}, 



dzUi 1 2=0 — OxUq \x=0, 

lim = 0. 
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Hence Ui is given by: 

We will show now that the profiles can be computed as any order. Assume 
that U*~ has been computed, is solution of the well-posed 1-D wave 
equation: 

' dttU^ - c^dxxUf = 0, (x, t) G]0, 7r[xM+, 

f/j^!x=o = u* 1 2=0- 
< f//U.=. = 0. 

C/+|i=o(3;) =0, {x>0}. 
^dtU+\t=Q = (), {x>0}. 
V^^i is then solution of the well-posed profile equation: 

' - c^d,,U*-^ = -dttU*-, {z < 0}, 

lim = 0. 

Let us answer the question asked in [11] : is bound to present boundary 
layer behavior in {x = 0^}, indeed its approximate solution is composed 
exclusively of boundary layer profiles, which describes quick transitions at 
the boundary using a fast scale in e. As a result of the loss in convergence 
induced by the boundary layer, the following estimate holds: 

\\U^ - f^||L2(]-oo,7r[xR+) = C'(e2). 

In [11], their small parameter is ;U = e^, as a result, adopting the same 
notations as them, our estimate writes: 

which is in agreement with the estimates given in [11] . Like in the penaliza- 
tion approach proposed by Bardos and Ranch p] and underlined by Droniou 
in [3] , the boundary layer only forms on one side of the boundary. The ap- 
proximation U^'^ of [/, is computed by taking [7'^+|^=o = U^~\x=o-, thus, in 
numerical applications, the boundary layer phenomenon also affects the rate 
of convergence of U^^ towards U, 
as e ^ 0"*". 
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